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ABSTRACT 

We explore the impact of spatial fluctuations in the intergalactic medium temperature on the Lya 
forest flux power spectrum near z ~ 3. We develop a semianalytic model to examine temperature 
fluctuations resulting from inhomogeneous H I and incomplete He II reionizations. Detection of these 
fluctuations might provide insight into the reionization histories of hydrogen and helium. Furthermore, 
these fluctuations, neglected in previous analyses, could bias constraints on cosmological parameters 
from the Lya forest. We find that the temperature fluctuations resulting from inhomogeneous H I 
reionization are likely to be very small, with an rms amplitude of < 5%, aTo/{To) < 0.05. More 
important are the temperature fluctuations that arise from incomplete He II reionization, which might 
plausibly be as large as 50%, crTb/(2o} ~ 0.5. In practice, however, these temperature fluctuations 
have only a small effect on flux power spectrum predictions. The smallness of the effect is possibly 
due to density fluctuations dominating over temperature fluctuations on the scales probed by current 
measurements. On the largest scales currently probed, k ~ 0.001 s km~^ (~0.1 h Mpc~^), the effect 
on the flux power spectrum may be as large as ^ 10% in extreme models. The effect is larger on 
small scales, up to ~ 20% at fc = 0.1 s km~^, due to thermal broadening. Our results suggest that 
the omission of temperature fluctuations effects from previous analyses does not significantly bias 
constraints on cosmological parameters. 

Subject headings: cosmology: theory - intergalactic medium - large scale structure of universe; quasars 
- absorption lines 



1. INTRODUCTION 

In the current theoretical picture of the Lya forest, 
most of the structure in the forest is a product of grav- 
itational instability. The absorbing gas is assumed to 
be in photoionization equilibrium with a spatially ho- 
mogeneous radiation field. On large scales the hydro- 
gen gas distribution follows the dark matter distribution, 
and on small scales it is Jeans p ressur e-smoothed (see 
ej?., "Cen et al."1994; 'Zhang et al."1995; 'Hernauist et al.' 
1996; Miralda-Escude et al. 1996; Mucckct et al. 1996; 
Bi & Davidson 1997; Bond & Wadslev 1997; Hui et al. 
1^ .Croft et al. .1998^ .Brvan et a,1.. .199ft .Dave et nh 
19991 iTheuns et all 119991: INusser & Haehneltl I1999D . 
This gravitational instability model of the forest, moti- 
vated by numerical simulations, seems to agree w ell with 
observa tions (e.g. [Cro ft ct al. 2002; McDonalc^_et_^ 
|2005bl: iTvtler et ain2004: Viel et al.. .2004 .Lidz et al. 

In this model, each Lya forest spectrum provides 
a one-dimensional map of the density field in the in- 
tergalactic medium (IGM). The Lya forest can thus 
be used to constrain the amplitude and slope of 
the linear matter power spectrum at z ^ 3 on 
scales of k ^ 0.1 — 5/iMp c~^ (see e .g.. iCroft et al.l 
[1998: McDonald et alJ l2f)f)fi lz^,lda ,rria, ga, et a,lJ 12001? 
ICroft et all 120021: iZaldarriaga et al.ir2003l) . When com- 
bined with measurements from Cosmic Microwave Back- 
ground (CMB) experiments and galaxy surveys, the Lya 
forest provides important constraints on cosmological pa- 
rameters (e.g. Scliak ct al. 2005; Viel ct al. 2005b). Re- 
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cently, very tight constraints on cosmological parame- 
ters were derived using measurements of the Lya forest 
flux power s pectrum from th e Slo an Digital Sky Sur- 
vey ( SPSS) llMcDonald eni\\2(M I2005H ISeliak et alJ 
|2005[) . The data samples used in the SDSS flux power 
spectrum measurements are almost two orders of mag- 
nitude larger than those used in previous measurements. 
The increase in statistical precision sets a high bar for 
the required control over systematic effects in theoretical 
predictions. The key issues are now to devise consistency 
checks for the gravitational instability model of the Lya 
forest, to quantify the accuracy of our theoretical model- 
ing, and to improve the modeling when possible. These 
steps are essential in order to estimate the systematic- 
error budget in the modeling, and to utilize the full sta- 
tistical power of the SDSS measurements. 

Toward this end, we point out that previous mod- 
els of the Lya forest adopt an over-simplified descrip- 
tion of the IGM reionization history and the result- 
ing thermal evolution. Specifically, previous models 
assume that reionization is sudden and uniform, so 
that in effect each gas element in the IGM experi- 
ences the same reionization history. In reality, reion- 
ization is likely to be an extended a n d inhomogeneous 
process llSokasian et al"l 120031 120041: iBarkana fc Loebl 
I20M iFurlanetto et al.ll2004l:lBabich &: Loehll200^ with 
some gas elements reionizing earlier than others, and 
hence cooling to lo wer temperatures by z ~ 3 (e.g. 
iHui fc Haimanll2003^1 . Furthermore, He II may be reion- 
ized by bright quasars, and the process may be incom- 
plete near z ~ 3. In this case, the IGM at z ~ 3 resembles 
a two-phase medium. The first phase consists of regions 
that have already been engulfed by the He III ionization 
fronts that are expanding around bright quasars. These 
regions, recently photo-heated by a hard quasar spec- 
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trum, may have te mperatures in excess of ~ 3 x IC* K 
l)Ahel fc Haehnelti ri999'). The second phase consists of 
regions where He II has yet to reionize, but have H I/He I 
reionized at early times. In this phase, gas elements with 
density near the cosmic mean will be significantly cooler, 
with temperatures of 1 x 10'* K. These temperature 
fluctuations should be imprinted in the Lya forest since 
the widths of Lya absorption lines, as well as the hy- 
drogen recombination coefficient, and hence the optical 
depth to Lya absorption, depend on temperature. 

The goal of this paper is to estimate the amplitude 
and spatial scale of the temperature fluctuations result- 
ing from H I and He II reionization, and to examine 
their impact on the flux power spectrum. Temperature 
fluctuations are particularly interesting because if they 
significantly impact the flux power spectrum, then their 
detection would likely provide insights into the reioniza- 
tion histories of hydrogen and helium. Furthermore, we 
reiterate that it is important to check whether omitting 
temperature fluctuations in the analyses will significantly 
bias the determination of cosmological parameters from 
the Lya forest. 

The outline of this paper is as follows. We begin with 
a brief overview of the theoretical model describing the 
Lya forest in § |21 We then estimate the amplitude of 
temperature fluctuations expected from inhomogeneous 
H I reionization in § |21 In § 0] we estimate the ampli- 
tude and scale of temperature fluctuations for a range of 
models describing He II reionization by bright quasars, 
using the observed quasar luminosity function as input. 
In § 121 we examine the impact of these fluctuations on 
the flux power spectrum. We conclude in §|Hland discuss 
possible future research directions. 

2. MODELING THE Lya FOREST 

In this section, we briefly review the standard theoret- 
ical model of the Lya forest in order to introduce no- 
tation, and highlight the approximations that we subse- 
quently test in t his paper. For m ore details, the reader 
can refer to e.g..' Hui et a 

The gas responsible for the absorption in the Lya 
forest is thought to be in photoionization equilibrium 
with a radiation background produced by star-forming 
galaxies and/or quasars. In this case, the abundance 
of neutral hydrogen scales like nni oc a(T)A^/r. Here, 
a(T) oc T^*''^ is the temperature dependent hydrogen re- 
combination coefficient, A is the baryon density in units 
of the cosmic mean, and F is the hydrogen photoioniza- 
tion rate. The optical depth to Lya absorption is propor- 
tional to the neutral hydrogen abundance (besides ther- 
mal broadening and peculiar velocities), which implies a 
simple power-law relationship between the Lya optical 
depth and the gas density. Therefore, if the gas temper- 
ature and F are known or can be modeled, then there 
exists a direct connection between the Lya absorption 
and the underlying density fluctuations. 

Indeed, the physics that sets the temperature of the 
absorbing gas is expected to be relatively simple. The 
gas temperature is determined largely by the competition 
between photoionization heating and adiabatic cooling 
(piralda-Escude & Recs 1994; Hui & Gncdin 1997). In 
this case, the temperature of the low density gas, where 
shock-heating should be unimport ant, is expected to b e 
tightly correlated with its density l)Hui &: GnediiJ[T99^ . 



In fact, these authors show that the gas temperature 
should be a power-law in the gas density: T — TqA"^^^. 
The numerical values of the power law index, 7, and 
the temperature at mean density, Tq, depend on when 
the gas wa s reioniz ed and the nature of the ion izing 
sources ('e.g.lHui fc Ctnedin 1997tlAbel Haehneltiil999.: 
Sokasian et al.l2002HTheuns et all2002aHHui fc Haimanl 

The standard assumption is that F, Tq and 7 are all 
spatially uniform, i.e. they have a single value through- 
out the entire IGM. To the extent that this is true, 
the neutral hydrogen density, and hence the Lya op- 
tical depth, scales as r = ArA-^~°-''^''~^\ where the 
proportionality constant Ar is independent of spatial 
position. In this simple case, fluctuations in the Lya 
forest absorption directly traces the underlying den- 
sity fluctuations. Several authors have investigated 
the possible 'contamination' from spatial fluctuations 
in the hydrogen photoionization rate, F, finding that 
these fluctuations should be quite small near z ~ 3 
llCroft et a].iri999t IMeiksin Whitel 120041 IcT^ l20(il 
iMcDonald et al.ll2005af) . Little attention, however, has 
been given to the assumption that Tq and 7 are also spa- 
tially uniform. 

In the case where Tq and 7 fluctuate spatially, the re- 
lation for the Lya optical depth generalizes to: 

T-A,A2-o-^(^-i)(l + J^J-o-^ (1) 

The new ingredient in this equation is the term 1 -I- Stq = 
Tq/ (Tq) which represents spatial fluctuations in the tem- 
perature of the gas at the cosmic mean density. Further- 
more, we consider spatial fluctuations in 7, i.e. 7 in this 
equation now depends on position. Additionally, fluctu- 
ations in the IGM temperature lead to spatial variations 
in the thermal broadening kernel, which will affect the 
Lya forest on small scales. 

The flux transmitted through the Lya forest, after tak- 
ing into account peculiar velocities and thermal broad- 
ening, is given by = e^'^. Fluctuations in the trans- 
mission are given by Sp = [F — {F))/(F), and the power 
spectrum oi Sp is the flux power spectrum, which we de- 
note by Pp (k) . Our goal then is to explore the impact of 
the temperature fluctuations encoded in Eq.^ and in the 
thermal broadening kernel, on the flux power spectrum 
Ppik). 

3. TEMPERATURE FLUCTUATIONS FROM 
INHOMOGENEOUS H I REIONIZATION 

We begin by considering the amplitude of tem- 
perature fluctuations resulting from H I reioniza- 
tion, which recent theoretical work has emphasi zed to 
be likely inhomogen eou s and e xtended (S okasian gLalJ 
2003; Barkana & Loe^ 120041: iFurlanetto et all 120041: 
iBabich fc Loeh .2006i) . We use the model of 
l^irfanetto'et'^Ll l)2004|) to describe the duration of ex- 
tende d H I reioniza tion, and scaling relations from 
Hui fc Gnedinl l)1997|) to estimate the amplitude of the 
resulting temperature fluctuations. Throughout this pa- 
per, we will work generally within the context of a model 
in which H I/He I are reionized at high redshift by star- 
forming galaxies, while He II is reionized close to z 3 
by quasars. This is not the only possibility (see Lidz et 
al. 2005, in prep., for a discussion), and so it is useful 
for us to explicitly separate out the temperature fluctua- 
tions that arise from H I reionization and those that arise 



Lya Forest Power Spectrum 



3 



N 

St 

T3 



0.4 
0.3 

0,2 

0.1 
0.0 




10 



15 



20 



Q. 
■D 



0.0100 



0.0010 



0.0001 




1.1 1.2 



1.3 



1.4 1.5 



1.6 



To [10' K] 



Fig. 1. — Panel (a): H I reionization redshift distribution for f 
= 6, 12, and 18. Panel (b): Temperature distribution at 2 = 3 for 
the same reionization models in (a). The thick (thin) curves are 
calculated using = 2.5 x 10* K (3 x 10* K). 



from He II reionization. The case that H I and He II are 
both reionized at high redshift will then be similar to our 
H I calculation: the only difference being that the tem- 
perature at reionization would likely be higher, and the 
duratio n of the reionizat ion process might be different. 

The iFurlanetto et alJ (|2004) model describes the 
growth and overlap of H II regions during reionization. 
The basic picture in this model is that large scale over- 
dense regions contain more ionizing sources, and are 
reionized earlier, than underdense regions. The model 
assumes that a galaxy of mass Wgai can ionize a mass cor- 
responding to ^rrigai, where C, is an unknown parameter 
describing how efficiently a g^alaxy can ionize surround- 
ing g as (see IFurlanetto fc Oh " 2005 and IFurlanetto et al.l 
l200fl for extensions to this model). A region is con- 
sidered ionized when the the fraction of mass in halos 
more massive than some minimum mass, mmin, exceeds 
a threshold set by the ionization efficiency of the sources: 
/coll > C~^- Here, fcow denotes the fraction of mass in the 
region which has collapsed into halos of mass larger than 
mmin. In this case, the size distribution of H II regions 
can be calculated in a similar manner to th e halo mass 
function in the excursion set formalism fe.g.. lBond et al.l 
ITOQlt lLacev fc Cole'lQQ.t). 

The Furlanetto et al. (2QQ^ model predicts, given 
the filling factor of ionized regions. It is given by 

Q = C/coU (2) 

where /coii is the global collapse fraction (different from 
/coll, the collapse fraction of a region with a given over- 
density). In the extended Press-Sch echter theory, the 
global collaps e fraction is given by ijBond et alJ Il99ll: 
ILacev fc Coy ri993): 



/coll = erfc 



V2<i„(^) 



(3) 



where 5c{z) is the critical density for collapse and (j'^^^{z) 
is the density variance on the scale of mmin, the minimum 
mass of an ionizing source. We take the minimum mass 
to be the mass corresponding to a virial temperature of 
10^ K, where atomic hydrogen line cooling is efhcient. 

The probability distribution of reionization redshift is 
related to the filling factor Q by 

dP _ dQ 

dzr dz 



(4) 



In Fig. HJi, we plot dP/dzr for C = 6, 12, and 18. The 
curves terminate when Q — 1, corresponding to the end 
of reionization. The figure clearly illustrates that H I 
reionization should be quite extended, with the whole 
process taking place over a Az of several. The end, and 
somewhat the duration, of reionization naturally depends 
on how efficiently the sources produce ionizing photons, 
which is quite uncertain. We therefore consider ^ = 6, 
12, and 18, in which case the end of reionization oc- 
curs at z = 6, 8, and 9 respectively. For our purposes, 
these choices of C bracket the interesting range of possi- 
bilities. It is uninteresting to consider sources that are 
much less efficient, beca use we k now the IGM is highly 
ionized below z ~ 6 (e.e. iFan et a l. 2002). On the other 
hand, more efficient sources would lead to a very early 
end to reionization, in which case the temperature at 
lower redshifts is completely insensitive to when precisely 
the gas was reionized owing to efficient Compton cooling 
|Hui & Haiman 2003). 

In order to investigate the temperature fluctuations 
that result from extended H I reioniz ation, we rely on 
analytic approximations bv lHui fc Gne din (1997) to de- 
scribe the thermal history of the IGM. These analytic 
approximations derive from the fact that the tempera- 
ture of the low density gas is primarily determined by 
photoionization heating and adia batic cooling . Under 
these simple physical conditions. iHui fc Gned in (199'3) 
give formulae for the evolution of Tq and 7 (see their 
Eq. 19 and 22). The input to these formulae is simply 
the temperature, T^, that a gas element reaches when it 
is reionized at redshift z^. Given T^, these formulae give 
the values of Tq and 7 at lower redshifts, z < Zr- 

The temperature at reionization, T^, is quite uncer- 
tain. It is determined both by the intrinsic spectrum 
of the ionizing sources, and the hardening of the spec- 
trum owing to absorption in the IGM ( Abel & Haehne^ 
1999). Scatter in the intrinsic and re-processed spectra 
will likely give rise to a distribution of T^. It is also 
unclear that gas elements reionized at different times, al- 
beit by sources with the same intrinsic spectrum, will 
reach the same temperature following reionization. We 
currently ignore these subtleties in our model, and as- 
sume that all gas elements will reionize to the same T^. 
The value = 2.5 x 10^ K is rea sonable assuming tha t 
galaxies are the ionizing sources IIHui fc Gnedinlll99l . 
We also investigate the effects of using a more extreme 
= 3 X lO'* K. Note that because of the high abun- 
dance of galaxies, each gas element will likely see the 
combined radiation from numerous sources during reion- 
ization. This will tend to average out the scatter in the 
intrinsic and re-processed spectra, so a uniform Tr is 
probably a good approximation. The same may not be 
true if quasars are the ionizing sources, since quasars are 
sparse and each gas element will only see the radiation 
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from one, or a few, nearby sources. 

With z and Tr fixed. To is a function of Zr only and we 
can find the distribution of To with the simple transfor- 
mation: 



dP 



dP 

dZr 



dTo 



dZr 



(5) 



The temperature distribution is plotted in Fig. ^p. The 
cutoff in the distribution at high To occurs because there 
is a definite end to reionization in our model, after which 
the probability of reionization is formally zero. The up- 
per temperature limit to the distribution is therefore set 
by the To of the most recently reionized gas elements. 
There is a sharp rise in the temperature distribution at 
low To, because gas elements reionized above z ~ 10 will 
have reached almost the same temperature at low red- 
shift, owing to efficien t Compton cooling at high redshift 
lIHui fc GnedinllT997l) . 

The temperature distribution at z = 3 spans a very 
narrow range in To, generally less than a few hundred 
degrees Kelvin, in all the models we consider. This is 
because the interplay between photoionization heating 
and adiabatic cooling drives the gas towards a thermal 
asymptote (Hui & Haiman 2003). Therefore, gas ele- 
ments that are reionized sufficiently early will approach 
the same asymptotic temperature at low redshifts. As a 
result, the temperature fluctuations at z = 3 owing to 
H I reionization are small. In the C = 6 model, with 
Tr = 2.5 X 10'' K, the level of temperature fluctuations is 
about 4%, <7T„/{To) = 0.04. In the C = 12 and 18 models, 
the levels of temperature fluctuations are even smaller at 
< 1%. Taking T^ = 3 x 10^ K does not give significantly 
larger temperature fluctuations. The amplitude of tem- 
perature fluctuations might be larger at higher redshifts, 
because the gas elements have less time to cool. How- 
ever, efficient cooling quickly causes the gas to approach 
the thermal asymptote, and we find that the amplitude 
of temperature fluctuations at z = 4 is very similar to 
that at z = 3. 

It is possible to have temperature fluctuations larger 
than what we have estimated. For instance, if He II 
is reionized alongside H I/He I, then temperatures of 
4 X 10"^ K or higher are plausible following reioniza- 
tion ijAbel fc Haehnel^ fl999). leading to larger temper- 
ature fluctuations. Furthermore, the analytic approxi- 
mation we employ assumes a constant spectrum for the 
ionizing background. In reality, the spectrum experi- 
enced by a gas element can be hardened significantly 
during reionization rela tive to the late time spectrum 
ijAbel fc Haeh nemil999ll . Therefore, for a given Tr, the 
analytic approximation tends to overestimate the late 
time temperature, since the hardened spectrum at reion- 
ization is assumed throughout. Processes such as recom- 
bination cooling are also neglected in the analytic ap- 
proximation, further contributing to the overestimation 
of the late time temperature. A more detailed calculation 
taking into account different cooling mechanisms and 
variations in the spectrum of the ionizing background will 
in general yield larger temperature fluctuations. Note 
also that the temperature fluctuations at z ~ 6 might be 
significantly larger than those at z ~ 3. This might bias 
constraints derived from z ^ 6 quasar spectra, such as 
the evolution o f the hydrogen photoionization rate (e.g. 
iFan et al.ll20(il . 



In summary, we find that while H I reionization can 
be quite extended, gas cooling erases temperature fluc- 
tuations over time. Therefore, temperature fluctuations 
from inhomogeneous H I reionization are likely negligi- 
ble at z = 3, especially if H I is reionized by z > 8. 
Even though possibilities such as a high temperature at 
reionization or evolution in the spectrum of the ioniz- 
ing background might lead to larger temperature fluctu- 
ations, it seems likely that there must be 'reionization 
activity' very near z ~ 3 in order for the temperature to 
fluctuate significantly at this redshift. 

4. TEMPERATURE FLUCTUATIONS FROM INCOMPLETE 
HE II REIONIZATION 

One possible scenario that can give rise to potentially 
large temperature fluctuations is when He II is reionized 
gradually by bright quasars. If He II reionization is still 
underway at z ~ 3, there should be hot He III bub- 
bles around bright quasars embedded in a much cooler 
background IGM, in which only H I/He I has reionized. 
The temperature in the hot bubbles may be as high 
as To ^ 3 X lO'^ K l)Abel fc Haehneltlll999|) . while the 
temperature outside of these bubbles should be around 
To ^ 10* K, as described in the previous section. In 
this case, the temperature fluctuations can be signifi- 
cantly larger than those arising from H I reionization. 
Additionally, the regions that recently reionized He II 
will be close to isothermal (7 — 1), while the cool ex- 
terior will have a steeper temperature-density relation, 
7 1.4 — 1.6, with the precise value depending on when 
H I/Hc I reionizcs. In this section, we begin by describ- 
ing a simple model for the growth of He III bubbles, and 
their subsequent thermal evolution. We will then discuss 
the resulting temperature distribution and power spec- 
trum of temperature fluctuations. 

4.1. Numerical Model 

Our procedure is to populate a simulation box with 
quasars by drawing sources from the observed quasar 
luminosity function (QLF). We then follow the growth 
of He III bubbles around each quasar source for a fixed 
quasar lifetime, tq, and record the subsequent thermal 
evolution inside the ionized regions. We can then mea- 
sure the statistics of the resulting temperature field, and 
use the temperature field, along with the density and pe- 
culiar velocity fields from a cosmological simulation, to 
study the statistics of the Lya forest. 

First, we will consider the time evolution of the He III 
ionized volume around an isolated quasar source. The 
most interesting period for temperature fluctuations is 
before the overlap of He III ionized regions is complete. 
In this pre-overlap phase, the growth of He III regions 
is simple to describe, since we can make the approxima- 
tion that a gas element will only see the radiation from 
the central quasar. The growth of the ionized region 
around the quasar is th en given by ijShapiro fc Girou3 
[19871 iMadaii et al.llT999|) : 



dV_ 
~dt 



N 



V 



tr. 



(6) 



In the above equation, V is the comoving volume of the 
ionized bubble, nfjo is the cosmic mean comoving num- 
ber density of helium atoms, and N is the number of 
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TABLE 1 

Model Parameters 



Model 


a(z = 2.1) 


a{z > 3.6) 


/3 


A 


B 


tq [yrs] 


Tr [K] 


Fiducial 


-3.28 


-2.58 


-1.78 


-7.31 


0.47 


10^ 


3 X 10'' 


Small Fluctuations (SF) 


-3.48 


-2.81 


-1.78 


-7.12 


0.32 


10^ 


3 X 10* 


Large Fluctuations (LF) 


-3.28 


-2.58 


-1.78 


-7.69 


0.77 


10^ 


4 X 10* 


Small Scale (SS) 


-3.28 


-2.58 


-1.78 


-7.69 


0.77 


10'5 


4 X 10* 


Large Scale (LS) 


-3.28 


-2.58 


-1.78 


-7.50 


0.77 


10* 


4 X 10* 


Large Flue, z = 4 (LF4) 


-3.28 


-2.58 


-1.78 


-7.12 


0.47 


10^ 


4 X 10* 


Note. — Parameters for the QLF 


(a, 0, A, 


and B) 


are based on 


the best fit values 



and uncertainties in IRichards et all i2005D and IFan et alJ i200lD . The normalization 
of the QLF from IRichards et alJ ll2005l) . = 5.96 x 10"''Mpc"^mag"\ is used in all 
models. 



He II ionizing photons emitted by the quasar per unit 
time. This equation assumes that all of the helium in 
the pre-He II reionized gas is singly ionized. Further, it 
assumes that all photons emitted above the He II thresh- 
old contribute to He II reionization, since absorption of 
He II ionizing photons by H I and He I is negligible ow- 
ing to their high ionization levels. The ionized region is 
assumed to be spherical, and grows according to Eq. 
while the quasar is active. The bubble is assumed to 
remain fixed at its final size after the quasar turns off, 
and the subsequent redshift evolution of To and 7 inside 
the bu bble is tracked by the formulae of iHui &: GnedinI 

In Eq. El a single volume averaged recombination rate 
is assumed for He III inside the bujable. The recombina- 
tion time is tree = ("-c.pCkHgjjjC) , where a^^^m is the 
recombination c oefficie nt to the excited states of He III 
(case B) fsee iMadau et a l. 1999), and fio.p is the proper 
electron number density. The clumping factor of He III 
is defined as C = ('^Hciii)/^Hcni- "^he value of C is 
quite uncertain, and values between C = 1 — 30 are com- 
monly used in the literature (see e.g. . IMadau et al1ll999t 
iMeiksin 2005). In our calculation we chose C = 1. Re- 
gardless, the effect of recombinations on the growth of 
He III bubbles should be limited, since the recombina- 
tion time tiec is much longer than the expected quasar 
lifetime tq. Assuming an IGM temperature of 3 x 10'' K, 
iioc is on the order of 10^ yrs at z = 3, much longer than 
tq, which we vary between 10^ — 10* yrs. A quasar of 
luminosity L will then be surrounded by a He III region 
with an approximate volume of 1^ oc Ltq at the end of 
its lifetime. 

The next ingredient in our modeling is our description 
of the abundance, spectrum, and lifetime of quasars. We 
para meterize the QLF with the standard double p ower 
law l|Bovle et al.llT988l IPeiHToM ICroom et alJl200l : 

We use measuremen ts of the QLF from 2SLAQ 
l|Richa,rds et a,].l l200l at 0.4 < z < 2.1, and SDSS 
l|Fan et al.ll200lh at z > 3.6. Note that the bright-end 
slope a from SDSS appears to be significantly different 
(> 2a) than that measured by 2SLAQ at low redshift. 
Also, the SDSS measurements at z > 3.6 do not probe 
the faint end of the QLF. In order to bridge the gap 
between the SDSS and 2SLAQ measurements, and to 



extrapolate to all relevant luminosities, we make several 
assumptions. First, at intermediate redshifts, we linearly 
interpolate between the low and high redshift bright-end 
slopes. Second, we use the faint-end slope, (3, and the 
normalization of the QLF, from 2SLAQ and assume 
that they remain fixed with redshift. Finally, we fixed 
by requiring that our model QLF matches the SDSS best- 
fit abundance of bright quasars. Specific ally, we match 
to the fitting formula of IFan et al.1 1)200 1|) ^. in which the 
number density of quasars with magnitudes Mi45o < — 26 
is given by log$(z, M1450 < -26) = A~ B{z - 3), with 
$ in units of Mpc~'^. 

Finally, we adopt a low luminosity cutoff for the quasar 
lumin osity function of imin — O.OlSi* fsee IMadau et alJ 
19991 for a d i scussi on), and the quasar spectrum from 
Madau et al.1 1)1999() . which goes as v~^-^ above the He II 
ionization threshold. All quasar sources are assumed to 
have the same spectrum and lifetime. 

With these considerations in mind, we investigate sev- 
eral models. Our strategy here is to span a conserva- 
tive range in the parameters that characterize the quasar 
sources, and hence a range in the amplitude and char- 
acteristic scale of the resulting temperature fluctuations. 
This is prudent given not only the observational uncer- 
tainties in these parameters, but also the uncertainties 
and approximations inherent in our modeling. The pa- 
rameters of these models are summarized in Table ^ In 
our fiducial model, we adopt the best fit values from 
2SLAQ and SDSS for the parameters in the QLF, and 
use tq = 10^ yrs and T,. = 3 x 10^ K for the quasar life- 
time and temperature at reionization, respectively. We 
then vary the parameters around these values, investi- 
gating models with small temperature fiuctuations (SF), 
and large temperature fiuctuations (LF), as detailed in 
Table In the LF model, in addition to varying the 
parameters of the QLF, we adopt a large temperature 
at reionization, T,. = 4 x 10'* K. In each case, we vary 
the parameters of the observed QLF within their allowed 
2-cr range. Finally, we vary the quasar lifetime in or- 
der to cover a range in the characteristic scale of the 
temperature fiuctuations. Specifically, the large scale 
(LS) temperature fluctuations model adopts tq — 10* 
yrs, while the small scale (SS) temperature fluctuations 

The cosmology adopted by Fan et al. ({Hm, f^Aj = {0.35, 
0.65, 0.65}) is slightly different from the one we use ({nm, ^K, 
= {0.3, 0.7, 0.7}), but we do not attempt to correct for this small 
difference. 
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model adopts tq = 10^ yrs. 

4.2. Temperature Distributions 

We construct realizations of the Tq and 7 fields in each 
of the models discussed in the previous section. We use 
a 80 Mpc comoving box, and 512'^ mesh points. This 
boxsize is convenient for overlaying on the cosmological 
simulation which we will describe in §|21 Throughout, we 
average the statistics in the fiducial, LF, and LS mod- 
els over five independent realizations, in order to reduce 
scatter owing to the small number of He III bubbles in 
our 80 Mpc simulation box. In the other models, the 
bubble distribution is sampled sufficiently well with a 
single realization. 

With the numerical models in hand, we proceed to 
measure the probability distribution of the temperature 
field. We first examine in Fig. [5^ the probability dis- 
tribution of He II reionization redshifts, in analogy with 
Fig. This figure is constructed by recording the red- 
shift at which each pixel in the simulation is first engulfed 
by an expanding He III bubble. We only plot our fidu- 
cial, SF, and LF models, since the other models have 
temperature distributions that are similar to that of the 
LF model. One can see that He II reionization is quite 
extended, and that these simple models are each con- 
sistent with incomplete He II reionization near z ~ 3. 
Our models a re essentially Monte-Carlo versions of the 
IMadau et al.l (^999) calculation and are consistent with 
these earlier calculations. 

The resulting temperature distributions at z = 3 are 
shown in Fig. \Bp. In our models, we assume that the 



background IGM, in which only H I/He I is reionized, 
has a uniform temperature at Tq = 10* K, as justified in 
§13 We additionally assume a uniform 7 = 1.4 for the 
background IGM. In the figure, one can clearly see the 
bimodal temperature distribution that results from in- 
complete He II reionization. The 7 distribution exhibits 
a similar bimodal distribution, with 7 ~ 1 inside He III 
regions, and 7 = 1.4 in the background IGM. The tem- 
perature distributions are quite broad, with r.m.s. fluc- 
tuation amplitudes of (TTo/(7o) = 0.24, 0.33, and 0.58 
(o-To = 4.8 X 10^ K, 6.4 X 10^ K, and 1.2 x lO'^ K) in 
the SF, fiducial, and LF models, respectively. Note that 
the higher temperatures of the hot regions in the LF 
model are not a consequence of that model's reionization 
history, but rather because we chose a larger T^ to max- 
imize fluctuations. Compared to the results from §13 we 
see that temperature fluctuations from incomplete He II 
reionization can be as much as a factor of 10 larger than 
that from extended H I reionization. 

The temperature distributions shown in Fig. |2h are 
each distributions at z = 3 in different models, but they 
also roughly represent the fluctuations during different 
stages of reionization. The filling factor of He III regions, 
Q, is indicated by the area under the hot component 
of the temperature distribution. In the LF model, Q is 
about 0.5 at z = 3. Here, the amplitude of temperature 
fluctuations is maximal, as the hot and cold regions oc- 
cupy comparable fractions of space. On the other hand, 
in the SF and fiducial models, reionization is consider- 
ably more complete. The temperature distribution is 
dominated by hot regions by z = 3, and the fluctua- 
tions are smaller. At a higher redshift. He II reionization 
is less complete in the fiducial and SF models, and the 
temperature distributions in these models will be similar 
to that in the LF model at z = 3. 

4.3. Power Spectrum of Temperature Fluctuations 

The temperature distributions presented above are in- 
formative, but they tell us nothing about the charac- 
teristic scale of the temperature fluctuations. To inves- 
tigate this, we measure the 3-dimensional power spec- 
trum of the temperature fleld, Tq, from our simulations. 
In Fig. we show the dimensionless power spectrum, 
A|,(fc) oc k^Prik), for each of the flve models we con- 
sider. Each power spectrum shows a well defined char- 
acteristic scale, with the power spectrum growing as k'^ 
on large scales, and falling as k~^ on small scales. The 
large and small scales trends in A'^{k) are a natural con- 
sequence of our model in which the source positions are 
uncorrelated (see below for comments on this approxi- 
mation). The assumptions that the bubbles have a well 
defined boundary, and that Tq within each bubble is uni- 
form, also contribute to the trends seen in A|,(fc). The 
characteristic bubble size is set by (see Eqs.|51and[7I) the 
break in the luminosity function, L,, the quasar spec- 
trum, and the quasar lifetime, tq. The characteristic 
comoving scale in the fiducial, SF, and LF models is 
R* ~ 5/i~iMpc, while it is i?* ll/i-^Mpc in the 
LS model, and i?* ^ 2h~^ Mpc in the SS model. 

Two simplifying assumptions in our model may cause 
the characteristic scale to be underestimated. First, our 
model does not take into account quasar clustering, and 
so we tend to underestimate the chance that ionized 
bubbles around neighboring quasars will overlap to form 
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large ionized regio ns. Assuming a q uasar bias bq — A 
(extrapolated from lCroom et al.l2005j) , we find that clus- 
tering enhances the average number of sources inside an 
ionized region by a factor of ~ 3 over that of a uniform 
distribution. However, even when the effect of clustering 
is included, there are on average < 1 additional active 
sources inside an ionized region of size i? ~ lQh~^ Mpc. 
Therefore, the clustering of active sources can safely be 
ignored. Quasar clustering will also lead to the clus- 
tering and overlap of fossil ionized regions. In this case, 
neglecting clustering may cause the typical volume of hot 
regions to be underestimated by as much as a factor of 3. 
It is therefore prudent to investigate models spanning a 
wide range of characteristic scales, as we have done. The 
second simplification in our model is in our treatment 
of multiple ionizing sources. When the ionization fronts 
of two or more bubbles overlap, the resulting large ion- 
ized region will expand so as to conserve ionizing photons 
coming from the multiple sources inside the combined re- 
gion. We neglect this subtlety in our modeling, noting 
that the probability for overlap is not significant around 
Q = 0.5, when temperature fluctuations are largest. 

We aim to explore the effects of temperature fluctua- 
tions on the statistics of the absorption in the Lya forest. 
Since in the standard picture of the forest most structure 
derives from gravitational instability, it is instructive to 
compare temperature fluctuations with density fluctua- 
tions. Which is a more important source of structure 
in the Lya forest, temperature or density fluctuations? 
We address this question in Fig. |31 where we include 
curves indicating the linear matter power spectrum, and 



the lPeacock fc Dodd s (1996) fit for the non-linear power 
spectrum. This may not be exactly the relevant com- 
parison: for instance, the optical depth to Lya absorp- 
tion scales with matter density as (1-1-5)^, while it only 
varies with temperature as (1 -I- 5to)~^''^ ■ Density fluctu- 
ations of a given amplitude are therefore amplified into 
larger optical depth fluctuations than temperature fluc- 
tuations of the same amplitude. Nonetheless, the com- 
parison is suggestive. Fig. O illustrates that on small 
scales, fc > 1 /iMpc~^, the density power is at least 1 — 2 
orders of magnitude larger than the temperature power 
in all models considered. On these scales, fluctuations 
in the density field will be the more important source of 
structure in the Lya forest. On the other hand, on large 
scales, the temperature power may be comparable, or 
even dominant, in comparison to the density power. For 
instance, in the LS model, the temperature power is ac- 
tually larger than the density power for k < 0.3 h Mpc~^. 

Our results seem to suggest that temperature fiuctu- 
ations may be an important source of structure in the 
Lya forest on large scales. There are, however, a few 
caveats to this intuition. First, as mentioned above, the 
optical depth scales more strongly with density than tem- 
perature. Second, current flux power spectrum measure- 
ments do not probe very large scales, as illustrated by 
the red dotted lines in Fig. O which show the range 
of scales probed by SDSS observations. This suggests 
that if the characteristic scale of the temperature fluc- 
tuations is large, much of the effect will be on scales 
larger than that probed by current measurements. Fi- 
nally, it is important to keep in mind, that the Lya for- 
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est provides a 1-d skewer through the IGM. As a result, 
the power spectrum of fluctuations along a line of sight, 
on large scales, includes aliased power from small wave- 
length modes transverse to the line of sight. This may 
tend to wash out some of the signal from temperature 
fluctuations on large scales. 

5. IMPACT OF TEMPERATURE FLUCTUATIONS ON THE 
FLUX POWER SPECTRUM 

We now arrive at the heart of our study: how do the 
temperature fluctuations, described in the previous sec- 
tion, impact the flux power spectrum? To answer this 
question, we turn to cosmological simulations. We com- 
bine the realizations of the Tq and 7 fields discussed in the 
previous section with the density and peculiar velocity 
fields from a cosmological simulation. We then extract 
absorption spectra from the simulation, and measure the 
flux power with and without temperature fluctuations. 

The cosmological simulation we use is a Hydro- 
Particle-Mesh (HPM) simulation (Gncdin & Hui"1998^ 
iHeitmann eTal. 20 05: Lidz ct al. 2006; Habib ct al. 2005, 
in prep.), with 2 x 5f2'' particles and 512^ mesh-points 
in an 80/i~^Mpc box. For our purposes, we want 
a simulation box that is large enough to sample the 
He III bubble distribution, while resolving the pressure- 
smoothing and thermal broadening scales. Our simu- 
lation represents a compromise between these require- 
ments. The resolution is inadequate for detailed predic- 
tions, as shown in the convergen ce studies presented in 
the Appendix of lLidz et al. (2006) . Furthermore, the ac- 
curacy of HPM has been criticized in the literature (e.g. 
iViel et al. 2005a). However, our present goal is only to 
investigate how the flux power spectrum differs with and 
without temperature fluctuations, rather than to make 
a detailed comparison with data. For this limited pur- 
pose, we believe our simulation is adequate. One fur- 
ther caveat is that the gas pressure force, as calculated 
in our HPM simulation, depends on the thermal history 
of the IGM. In principle, we should use the fluctuat- 
ing temperature field in our models to calculate the gas 
pressure in our HPM simulation. However, we ignore 
this effect and use HPM simulations calculated with a 
uniform temperature, including the temperature fluctu- 
ations only when we construct the artificial Lya forest 
spectra. The flux power spectrum, particularly on SDSS 
scales, depends rather w eakly on gas pressure-sm oothing 
llMcDonald et al.ll2005hl: IViel & Haehneltll2006D . and so 
this is probably a good approximation. 

We extract artificial spectra at z = 3 from the HPM 
simulation box in the usual way (see Eq. ^ , incorporat- 
ing peculiar velocities and thermal broadening. This is 
done for each of our five models, and we compare the flux 
power spectrum in the fluctuating temperature model 
with that in a similar model without temperature fluctu- 
ations. In the models without temperature fluctuations. 
To and 7 are each set to their global averages in the 
corresponding models that include temperature fluctu- 
ations. In each case, we normalize the quantity Aj- in 
Eq. n to match the observed mean transmitted flux at 
z — 3, which we take to be (F ) = 0.68, close to the value 
measured bvlMcDonald et all ((TOOO). 

The results of this calculation are shown in Fig. 0] and 
Fig. |31 where we parameterize the effect of temperature 
fluctuations by the fractional difference in the flux power 



spectrum between a fluctuating temperature model and 
a no fluctuations model. The first feature to notice is 
simply that the effect of temperature fluctuations on the 
flux power spectrum is quite small on all scales examined. 
Closer examination reveals that the flux power spectrum 
is boosted on large scales and on small scales compared 
to equivalent models without temperature fluctuations, 
while there is a very slight suppression on intermediate 
scales. We will discuss each effect in turn. 

The first effect is due to increased structure in the 
forest on large scales, contributed by the hot, isother- 
mal He III bubbles. The flux power is boosted on 
large scales because the flux transmission is sensitive 
to the temperature-dependent recombination coefficient, 
and the spatially fluctuating temperature-density rela- 
tion. In Fig. 01 concentrating on large scales for the 
moment, we study the effect for models with varying lev- 
els of temperature fluctuations: the flducial, SF, and LF 
models. Specifically, we compare the fractional difference 
between the flux power spectrum with and without tem- 
perature fluctuat ions and the 1-a error b ars on the SDSS 
measurements of lMcDonald et ahl l)2006j) . The fractional 
difference between the models is always smaller than the 
1-a SDSS errors. The difference is only at the 5% level 
on the largest scales probed in the LF model. The effect 
also appears to diminish rather quickly with decreasing 
temperature fluctuation strength, as illustrated by the 
other two models in the figure. 

In Fig. we show the same comparison for models 
in which we vary the characteristic scale of the temper- 
ature fluctuations: the SS, LF, and LS models. Each 
of these models has approximately the same fluctuation 
strength, and only differ in their characteristic scale. In 
this case, the effect can be as large as ~ 10% on the 
largest scales probed, comparable to the SDSS 1-a error 
bars at this redshift. In the models where the character- 
istic scale of the temperature fluctuations is smaller, the 
effect on the flux power spectrum is smaller. This ap- 
pears to be consistent with the interpretation suggested 
in § 14.31 density fluctuations generally swamp temper- 
ature fluctuations, unless the temperature fluctuations 
have a large characteristic scale. 

The next effect is a boost in the small scale power, 
which is the result of thermal broadening. The figures 
illustrate that the models with temperature fluctuations 
typically have ^ 20% more power on scales of k ^ 0.1 
s/km than corresponding models with a uniform To and 
7. In Fig. 01 and we compare this boost in small scale 
power t o the 1-a statistica l error-bars on the measure- 
ment of iCroft et alJ l)2002fl , which is the most precise 
measurement to date on these scales. The enhanced 
small scale power in the models with temperature fluc- 
tuations can be understood as follows. On small scales, 
the power spectrum in a fluctuating temperature model 
will approximately be a filling-factor weighted average of 
the power spectrum of hot regions and that of cold re- 
gions. Roughly speaking, the power spectrum on small 
scales is e xponentially suppressed with increasing tem- 
perature ijZaldarriaga et alJl2001a|) . The weighted av- 
erage we mention, and hence the power spectrum on 
small scales in fluctuating temperature models, is there- 
fore dominated by the cold regions. The fluctuating tem- 
perature model will then have more small scale power 
than a uniform temperature model with the same mean 
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Flux power spectrum with and without temperature fluctuations 
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temperature as the fluctuating model. 

Finally, there is a slight suppression in the flux power 
on intermediate scales, typicaUy around a few percent. 
This suppression occurs on scales in which the temper- 
ature field has little power, and on scales too large for 
thermal broadening to have an effect. The suppression 
is likely a consequence of the fact that the normalization 
Ar required to match the observed mean transmitted flux 
is slightly smaller in models with a fluctuating tempera- 
ture field. The smaller Ar in the fiuctuating temperature 
models implies that density fiuctuations of a given ampli- 
tude are translated into lesser optical depth fluctuations 
on intermediate scales. 

One might expect temperature fluctuations to have a 
larger effect at high redshift, when the amplitude of den- 
sity fluctuations is smaller. We test this by considering 
a model that, at z = 4, has very similar temperature 
fluctuations to those in our LF model at z = 3. The pa- 
rameters of this model, which we call LF4, are detailed 
in Tabled We flnd that the effect, again parameterized 
by the fractional difference between the models with and 
without temperature fluctuations, is only a couple per- 
cent larger on large scales in the LF4 model. Even though 
the effect is larger at z = 4, it is less significant in the 
sense that the statistical errors on the SDSS measure- 
ment are larger at z — A than at z = 3. 

6. CONCLUSION 

In this paper, we have estimated the level of tem- 
perature fluctuations expected in the IGM at z ^ 3 
from extended H I reionization and incomplete He II 
reionization. We find that the temperature fluctuations 
from extended H I reionization should be quite small, 
CTo/iTo) ^ 5%, while the fiuctuations from incomplete 



Fig. 5.— Top: Similar to Fig.glbut for the LF, SS, and LS mod- 
els. Bottom: Flux power spectrum with and without temperature 
fiuctuations in the LF model. 



He II reionization might be as large as ^ 50%. These 
fluctuations should have only a small effect on the flux 
power spectrum: on large scales, k ~ 0.001 s/km, tem- 
perature fluctuations lead to an increase in the 2 = 3 
flux power spectrum by at most ~ 10%. On small scales, 
fc ~ 0.1 s/km, fluctuations in the thermal broadening 
scale boost the power by ~ 20%. 

Further study is required to quantify the effects of 
temperature fluctuations on cosmological parameter con- 
straints. In particular, note that in the LS model (see 
Fig. EJ, the effect of temperature fluctuations are at, 
or larger than, the 1-cr level over several independent 
data points. A detailed investigation will require a multi- 
paranieter fit to the observed data (see iMcDonald et aO 
l2005bll . which is outside the scope of the present paper. 
However, we can anticipate the results of a more de- 
tailed investigation by comparing with the effects of UV 
background fiuctuations, as studied in ' McDonald et all 
(2005a b). Fluctuations in the UV background and tem- 
perature fluctuations have a similar influence on the am- 
plitude and shape of the flux power spe ctrum, at l east 
on large scales ( see Fig . 12 oflMcDonal d et al.l 12005b). 
IMcDonald et all lj2005bll included the effect of UV back- 
ground fluctuations in a multi-parameter fit, and found 
little effect (< 1%) on the inferred values of the ampli- 
tude and slope of the linear power spectrum. The general 
reason for this insensitivity is that the effective errors on 
the flux power spectrum are la,rger than the raw statis- 
tical errors on the data, after McDonald et al. (2005h|) 
marginalize over other, more significant, effects. Simi- 
larly, we do not expect temperature fluctuations to have 
a substantial effect on the values of the amplitude and 
slope of the linear power spectrum. 

For the purpose of cosmological parameter estimation, 
it would be instructive to know whether the effect of tem- 
perature fluctuations is degenerate with any parameters 
describing the Lya forest. As we mentioned in the pre- 
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vious paragraph, the effect of temperature fluctuations 
may be degenerate with those of UV background fluctua- 
tions on large scales. However, they have different effects 
on small scales. A d ditiona lly, inspect ing, e.g. Fig. 13 
of McDon ald e1~an l|2n05bD (see also iViel fc Haehneiil 
1200 6') . shows that the boost in the flux power on large 
and small scales expected from temperature fluctuations 
is not closely mimicked by changing any single model- 
ing parameter. In principle, this means that the effect 
of temperature fluctuations is likely distinguishable from 
other effects. The redshift evolution of the effect poten- 
tially provides an additional diagnostic. However, even 
though the number of quasar spectra will likely more 
than double by the time SDSS is finished, detecting tem- 
perature fluctuations in the flux power spectrum will re- 
main a challenge owing to the smallness of the effect. 

Numerous improvements could be made in our simple 
modeling. First, we place quasars at random positions in 
our simulation box, which ignores quasar clustering. This 
effect is already discussed in § 14.31 Here, we note that 
since quasars are very sparse sources, ignoring source 
clustering is a much better approximation than in the 
case that the sources are galaxies, i n which case t his ap- 
proximation is quite poor (Furlanet to et al.ll2004|) . Fur- 
thermore, in reality, quasars should reside in very mas- 
sive halos, rather than at random positions. As a result, 
we ignore the fact that very close to the quasar the gas 
will be overdense, tending to cancel out the enhanced 
transmission owing to the hot He III bubble around the 
quasar. The bubbles are, however, ^ 10 Mpc in size, 
over which the density contrast will be small. 

In modeling temperature fluctuations (the 'thermal 
proximity effect'), we also ignored the 'radiation prox- 
imity effect'. In reality the intensity of ionizing radia- 
tion will be enhanced over that of the radiation back- 
ground close t o an active quasar (e.g. IScott ct al. 2002; 
iRpflinde et al] r2Q05l. Nearby dead quasars, 'light cchos' 
of enhanc ed radiation will remain, propagating out into 
the IGM l)Croftll2004|) . These effects are ignored in our 
modeling, and the radiation background is treated as uni- 
form. In any case, the effects of the radiation proximity 
effect should be small compared to that from tempera- 
ture fluctuations. This is because the radiation proximity 
effect has a characteristic time scale of tq ~ 10^ yrs, short 
compared to the characteristic time scale of temperature 
fluctuations, ^cooi ~ 10^ yrs. Here, we define the cool- 
ing time tcooi to be the time it takes for a gas element 
to cool to half its original temperature assuming only 
adiabatic cooling. The longer characteristic time scale, 
and the fact that temperature fluctuations and the radi- 
ation proximity effect have similar physical characteristic 
scales, imply that a much larger fraction of space is af- 
fected by temperature fluctuations than by the radiation 
proximity effect. 

Additionally, we used a simplifled 'light bulb' model 
for quasar activity, in which each quasar shines at con- 
stant luminosity for the duration of its lifetime, tq, 
and emit no light thereafter. In reality, this is prob- 
ably a poor approximation to the quasar light curve, 
since quasars likely spend extended periods at low 
luminosity, going into or coming out of their peak 
luminosity phase ("Hopkin s et al.l T200!tI : iSpringel et al.l 
[2005.1. Finally, qua sars r aay launch large outflows (e.g. 
ISr.aTma,nieco fcOhll20fll IPi Matteo et a1.l 12001. which 



could also modify the absorption in their vicinity. 

In spite of all of these possible complications, we strove 
to cover a wide range in the amplitude and scale of tem- 
perature fluctuations in our analysis. It is unlikely that 
the problems discussed above will lead to an effect that 
lies far outside the range probed by our study. It should 
therefore be a fairly secure conclusion that temperature 
fluctuations do not significantly impact the Lya forest 
flux power spectrum on large scales. Our results provide 
additional support for the Lya forest as a robust probe 
of cosmology. 

Even though temperature fluctuations seem to have a 
small effect on the flux power spectrum, they may have 
a larger effect on other observables of the Lya forest. 
The flux power spectrum, while being the best measured 
statistic in the Lya forest, may be a poor statistic to use 
in searching for temperature fluctuations. There have 
been several searches for temperat ure fluctuations in the 
Lya forest with w avelet analyses l|Theuns et al.ll2002bl: 
lZaldarriagall2003) . These searches have failed to detect 
temperature fluctuations, but have been carried out us- 
ing only very small data samples. It would be interest- 
ing to investigate the expected signal from these searches 
given our modeling. A related statistic that might be sen- 
sitive to temperature fluctuations is the Lya forest 'tri- 
spec trum', as d efined in Zaldarriaga ct al. ( 2001b) (see 
also iFang fcW' hitc 2004). This statistic measures the 
scatter in the small scale power, as a function of scale, 
from region to region in the IGM. We might, therefore, 
expect a larger tri-spectrum in our models with temper- 
ature fluctuations than in models with no temperature 
fluctuations. 

One final point is that if He II reionization is under- 
way at 2 = 3, this might cause biases in estimates of 
the IGM temperature from the flux power spectrum, 
and the ionizing background derived from the Lya for- 
est proximity effect. The enhancement we find in the 
flux power spectrum on small scales might imply a slight 
bias in measurements of the IGM tem perature from the 
small scale flux power spectrum (e.g. iZaldarriaga et all 
2001a), and other similar measurements. A crude esti- 
mate of the bias is as follows. Approximately, thermal 
broadening suppresses the flux power exponentially so 
that Ppik) oc exp(— fc^cr^), where tr^ = khT/mH, and 
fcb is Boltzmann's constant. We found that temperature 
fluctuations increase Ppik) at fc = 0.1 s/km by 20%, 
which thereby implies a ^ 10% underestimate of the 
temperature in models that assume a uniform temper- 
ature. We caution, however, that the mean temperature 
is an incomplete description of the IGM thermal state in 
the presence of inhomogeneous reionization and/or large 
quantities of hot, shocked gas. Temperature fluctuations 
may also bias estimates of the ionizing background from 
the quasar proximity effect, which assumes that the IGM 
is at the cosmic mean temperature close to the quasar. 
We are investigating this, and other possible biases in 
the constraints from the proximity effect, using radiative 
transfer simulations. 

AL thanks Katrin Heitmann and Salman Habib for 
their collaboration in producing the HPM simulation 
used in this analysis. We thank the anonymous referee 
for useful comments on our manuscript. We also thank 
Scott Buries, Steve Furlanetto, John Huchra, and Peng 



Lya Forest Power Spectrum 



11 



Oh for useful discussions on these, and related topics. 
This work was supported in part by NSF grants ACI 
96-19019, AST 00-71019, AST 02-06299, AST 03-07690, 
and NASA ATP grants NAG5-12140, NAG5-13292, and 



NAGS- 13381. Some of the simulations were performed 
at the Center for Parallel Astrophysical Computing at 
the Harvard-Smithsonian Center for Astrophysics. 



REFERENCES 



Abel, T. & Hachnclt, M. G. 1999, ApJ, 520, L13 
Babich, D. & Locb, A. 2006, ApJ, 640, 1 
Barkana, R. & Locb, A. 2004, ApJ, 609, 474 
Bi, H. & Davidsen, A. F. 1997, ApJ, 479, 523 
Bond, J. R., Cole, S., Efstathiou, C, & Kaiser, N. 1991, ApJ, 379, 
440 

Bond, J. R. & Wadslcy, J. W. 1997, in ASP Conf. Ser. 

123: Computational Astrophysics; 12tli Kingston Meeting on 

Theoretical Astrophysics, 323 — h 
Boyle, B. J., Shanks, T., & Peterson, B. A. 1988, MNRAS, 235, 

935 

Bryan, G. L., Machacek, M., Anninos, P., & Norman, M. L. 1999, 
ApJ, 517, 13 

Cen, R., Miralda-Escude, J., Ostriker, J. P., & Ranch, M. 1994, 

ApJ, 437, L9 
Croft, R. A. C. 2004, ApJ, 610, 642 

Croft, R. A. C, Hu, W., & Dave, R. 1999, Physical Review Letters, 
83, 1092 

Croft, R. A. C, Weinberg, D. H., Bolte, M., Buries, S., Hernquist, 
L., Katz, N., Kirkman, D., & Tytler, D. 2002, ApJ, 581, 20 

Croft, R. A. C, Weinberg, D. H., Katz, N., & Hernquist, L. 1998, 
ApJ, 495, 44 

Croom, S. M., et al. 2005, MNRAS, 356, 415 

Croom, S. M., Smith, R. J., Boyle, B. J., Shanks, T., Miller, L., 
Outram, P. J., & Loaring, N. S. 2004, MNRAS, 349, 1397 

Dave, R., Hernquist, L., Katz, N., & Weinberg, D. H. 1999, ApJ, 
511, 521 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

Fan, X., Narayanan. V. K., Strauss, M. A., White, R. L., Becker, 

R. H., Pentericci, L., & Rix, H.-W. 2002, AJ, 123, 1247 
Fan, X., et al. 2001, AJ, 121, 54 
Fang, T. & White, M. 2004, ApJ, 606, L9 

Furlanetto, S. R., McQuinn, M., & Hernquist, L. 2006, MNRAS, 

365, 115 

Furlanetto, S. R. & Oh, S. P. 2005, MNRAS, 363, 1031 
Furlanetto, S. R., Zaldarriaga, M., & Hernquist, L. 2004, ApJ, 613, 
1 

Gnedin, N. Y. & Hui, L. 1998, MNRAS, 296, 44 

Heitmann, K., Ricker, P. M., Warren, M. S., & Habib, S. 2005, 

ApJS, 160, 28 

Hernquist, L., Katz, N., Weinberg, D. H., & Miralda-Escude, J. 

1996, ApJ, 457, L51+ 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B., & Springel, V. 2005, ApJ, 630, 716 
Hui, L. & Gnedin, N. Y. 1997, MNRAS, 292, 27 
Hui, L., Gnedin, N. Y., & Zhang, Y. 1997, ApJ, 486, 599 
Hui, L. & Haiman, Z. 2003, ApJ, 596, 9 
Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 
Lidz, A., Heitmann, K., Hui, L., Habib, S., Ranch, M., & Sargent, 

W. L. W. 2006, ApJ, 638, 27 



Madau, P., Haardt, F., & Rces, M. J. 1999, ApJ, 514, 648 
McDonald, P., Miralda-Escude, J., Ranch, M., Sargent, W. L. W., 

Barlow, T. A., Cen, R., & Ostriker, J. P. 2000, ApJ, 543, 1 
McDonald, P., Seljak, U., Cen, R., Bode, P., & Ostriker, J. P. 2005a, 

MNRAS, 360, 1471 
McDonald, P., et al. 2005b, ApJ, 635, 761 
McDonald, P., et al. 2006, ApJS, 163, 80 
Meiksin, A. 2005, MNRAS, 356, 596 
Meiksin, A. & White, M. 2004, MNRAS, 350, 1107 
Miralda-Escude, J., Cen, R., Ostriker, J. P., & Rauch, M. 1996, 

ApJ, 471, 582 

Miralda-Escude, J. & Rees, M. J. 1994, MNRAS, 266, 343 
Muecket, J. P., Petitjean, P., Kates, R. E., & Riediger, R. 1996, 

A&A, 308, 17 
Nusser, A. & Haehnelt, M. 1999, MNRAS, 303, 179 
Peacock, J. A. & Dodds, S. J. 1996, MNRAS, 280, L19 
Pei, Y. C. 1995, ApJ, 438, 623 
Richards, G. T., et al. 2005, MNRAS, 360, 839 
RoUinde, E., Srianand, R., Theuns, T., Petitjean, P., & Chand, H. 

2005, MNRAS, 610 
Scannapieco, E. & Oh, S. P. 2004, ApJ, 608, 62 
Scott, J., Bechtold, J., Morita, M., Dobrzycki, A., & Kulkarni, V. P. 

2002, ApJ, 571, 665 
Seljak, U., et al. 2005, Phys. Rev. D, 71, 103515 
Shapiro, P. R. & Giroux, M. L. 1987, ApJ, 321, L107 
Sokasian, A., Abel, T., & Hernquist, L. 2002, MNRAS, 332, 601 
Sokasian, A., Abel, T., Hernquist, L., & Springel, V. 2003, MNRAS, 

344, 607 

Sokasian, A., Yoshida, N., Abel, T., Hernquist, L., & Springel, V. 

2004, MNRAS, 350, 47 
Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 

776 

Theuns, T., Leonard, A., Schayc, J., & Efstathiou, G. 1999, 

MNRAS, 303, L58 
Theuns, T., Schayc, J., Zaroubi, S., Kim, T., Tzanavaris, P., & 

Carswell, B. 2002a, ApJ, 567, L103 
Theuns, T., Zaroubi, S., Kim, T., Tzanavaris, P., & Carswell, R. F. 

2002b, MNRAS, 332, 367 
Tytler, D., et al. 2004, ApJ, 617, 1 
Vicl, M. & Hachnclt, M. G. 2006, MNRAS, 365, 231 
Viel, M., Hachnclt, M. G., & Springel, V. 2004, MNRAS, 354, 684 
— . 2006, MNRAS, 367, 1655 

Viel, M., Lesgourgues, J., Haehnelt, M. G., Matarrese, S., & Riotto, 

A. 2005b, Phys. Rev. D, 71, 063534 
Zaldarriaga, M. 2002, ApJ, 564, 153 

Zaldarriaga, M., Hui, L., & Tegmark, M. 2001a, ApJ, 557, 519 
Zaldarriaga, M., Scoccimarro, R., & Hui, L. 2003, ApJ, 590, 1 
Zaldarriaga, M., Seljak, U., & Hui, L. 2001b, Ap.l, 551, 48 
Zhang, Y., Anninos, P., & Norman, M. L. 1995, ApJ, 453, L57-I- 



